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We have performed experimental quantum state tomography of NOON states with up to four 
photons. The measured states are generated by mixing a classical coherent state with spontaneous 
parametric down-conversion. We show that this method produces states which exhibit a high fidelity 
with ideal NOON states. The fidelity is limited by the overlap of the two-photon down-conversion 
state with any two photons originating from the coherent state, for which we introduce and measure 
a figure of merit. A second limitation on the fidelity set by the total setup transmission is discussed. 
We also apply the same tomography procedure for characterizing correlated photon hole states. 
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PACS numbers: 42.50.-p,42.50.Dv,42.50.Xa,42.50.Ar 



I. INTRODUCTION 

Quantum state tomography (QST) is a well known 
method for full characterization of experimentally formed 
quantum states. It is common in experimental quan- 
tum optics to generate multiphoton quantum states in 
which each spatial mode is occupied by exactly one pho- 
ton [J, [3] • More generally, one could consider multipho- 
ton states in which each o ptic al mode may accommodate 
more than one photon [3H12||. 'NOON states' which are 
N photon states in which all the photons occupy cither 
one mode or the other are an example of such states. The 
name 'NOON' stems from their mathematical Fock state 
notation [H [3] 

\NOON) A , B = -j=(\N, 0)a,b + |0, N) A , B ) (1) 

where A and B denote two optical modes. It was recently 
shown that NOON states can be generated with large 
photon numbers (high-NOON) [tJ, by mixing a classical 
coherent state (CS) and quantum spontaneous paramet- 
ric down-converted (SPDC) light on a standard beam- 
splitter [IH, [l6| . The resulting states were measured by 
showing super-resolution in the A-fold coincidence mea- 
surement. However, full characterization of the result- 
ing states requires QST. In this paper we experimen- 
tally implement a method for performing QST measure- 
ments on high photon number states in two modes. The 
method is based on a convex optimization algorithm for 
the maximum-likelihood function. In addition, we apply 
the method to reconstruct the density matrices for an- 
other class of states, the correlated photon holes (CPH) 
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IL MATHEMATICAL ANALYSIS 

The density matrix of a general mixed-state of N in- 
distinguishable photons in two modes, assumed here to 
be the two orthogonal polarization modes of one spa- 
tial mode, has (N + l) 2 elements that correspond to all 
-/Vth-order coherences. These Ath-order coherences are 
expressed by fl~7j 



(a\ n al N - n a?a"- m ) 



(2) 



where ai and a 2 are the two operators of the two orthog- 
onal modes of the electric field. It was recently shown 
that by measuring only A-fold coincidence, it is possible 
to determine all Nth. order coherences for arbitrary N, 
using SU(2) transformations [Hj]. Using Jones calculus 
notation, we can parameterize this transformation, using 
two angles (9 and <fi), as 



(3) 



where b\ and bi are the two orthogonal operators of the 
electric field, representing the state following the linear 
transformation U(0,<p). This transformation can be im- 
plemented in many equivalent ways 0, [l8[ , one of which 
includes using a liquid crystal (LC) as a retarder of phase 
cj), followed by a half wave plate (HWP) with an orienta- 
tion angle of 9 + 7r/4 with respect to the LC, 
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The Nth-order intensity moments of polarization mode 
Si, (b\ b^)k, abbreviated as Bk, can be measured while 
the LC and HWP are set to 4>k and 9k respectively, and 
the measurement is iterated for (N + l) 2 different angle 
pairs (meaning k = 1,2, (N + l) 2 ). Using Eqs. 
and ([3]) , we can relate the measurement of Bk to all the 
A^th-order coherences forming a set of (N + l) 2 linear 
equations: 



B k = R 



k Pn 



(5) 



where m and n are by convention summed from to N, 
andi?r = (-l) m -72 Ar C)Ox 



(cos ( 



sin6' fe )' l + m (cos6i /£ + sin0 fc ) 2JV - n - m eW m - n ). 



The linear set of equations in Eq. (J5) has a solution 
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FIG. 1: (Color online) Experimental setup for quantum 
state tomography of a multiphoton state in two modes, (a) 
Schematic of the setup depicting a prepared beam forming 
a multiphoton quantum state in two orthogonal polarization 
modes ai and cb, undergoing a linear state transformation 
introduced by a half wave plate (HWP) and a liquid crystal 
(LC). The transformed state is photon- number resolved, in 
one of its polarization modes b\. (b) Detailed layout of the 
setup. A pulsed Ti:Sapphire oscillator with 120 fs FWHM 
pulse width and 80 MHz repetition rate is doubled using a 
2.74mm LBO crystal to obtain 404nm ultra-violet pulses with 
maximum power of 225mW. These pulses then pump collinear 
degenerate type-I SPDC at a wavelength of 808nm using a 
1.78mm long BBO crystal. The SPDC (H polarization) is 
spatially and temporally overlapped with a CS (V polariza- 
tion) using a polarizing beam-splitter cube (PBS). A ther- 
mally induced drift in the relative phase, ip ca , between H and 
V polarizations is corrected every few minutes using a LC 
phase retarder. The coherent light amplitude is adjusted us- 
ing a variable attenuator. Modes &i ,&2 are realized collinearly 
at the ±45° X, Y polarizations respectively. These modes 
then pass through a LC aligned with the X, and Y polariza- 
tions introducing a relative phase shift <j) between them, and 
a HWP rotated by 9 in the H and V plane. Modes b\, hi 
are oriented at polarization axes H and V respectively. The 
spatial and spectral modes are matched using a PMF and a 
3nm FWHM band-pass (BPF) filter centered at 808nm. Pho- 
ton number resolving detection is performed using multiple 
single photon counting modules (Perkin Elmer). Additional 
components: Dichroic mirrors, short-pass filter (SPF). 



from which p mn can be uniquely determined when the 
equals (TV + l) 2 . This condition can be 
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TABLE I: Angles of HWP rotation and LC phase in QST 
measurements for a specific N. Each QST for a fixed number 
of photons, N, constitute a set of consecutive TV-fold coin- 
cidence measurements in polarization mode &i, where each 
measurement in this set is done while the HWP and LC are 
prepared in a pair combination of 8 and </> from this table. 
Each QST measurement set is composed out of all (N + l) 2 
pair combination in this table (except for odd TV). 
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"In the case of odd TV, a measurement for the pair 
= -%,<j> = should be added to all other (TV + l) 2 - 
pair combinations from the table. 



k = 1, (TV + l) 2 and for each such pair the TV- fold 
coincidence Bk is measured. An Example for possible 
sets are given for different TV in Table fl] 

Determination of a quantum state is implemented by 
reconstruction of the density matrix from a set of mea- 
surements and the set of linear equations in Eq. ([5]) re- 
lating the measurements to the coherences. A naive way 
to perform this reconstruction is by forming p and B as 
vectors, and R (TV + l) 2 x (TV + l) 2 matrix with 
columns indexed by mn and rows by fc and inverting 
Eq. (J5]) to find the density matrix: p mn = (i?™") _1 i?fe. 
This can be problematic as the resulting density matrix 
can be non-physical due to experimental noise, for exam- 
ple the estimated density matrix may well not generally 
be semi-positive definite. The most probable physical 
density matrix is obtained by employing the method of 
maximum-likelihood; this estimation technique is a well 
known method in statistics that is used in QST to deter- 
mine the most probable density matrix compliant with a 
set of measurements . 

£(p e st), denned as 



£(Pest) = 



(N+l) 

E 

k=l 
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met by choosing a set of (TV + 1) pairs 



c ) , where 



is a convex probability function that an estimated density 
matrix p est would produce the set of measurements of N- 
fold coincidence, Bk- 

This function was checked using a Monte-Carlo (MC) 
method to be a consistent, unbiased and efficient likeli- 
hood function estimator of the density matrix (for dif- 
ferent input states, with and without simulated noise). 
The problem of quantum state estimation for a given N 
is now reduced to finding the minimum for the function 
£-(Pest) using an optimization algorithm. It was noted be- 
fore [2fJ| that the likelihood function is convex and hence 
its global minimum can be found by a convex optimiza- 
tion algorithm. To perform convex optimization we used 
the free optimization package YALMIP [2l| under the 
physical constraints on the density matrix of unit trace, 
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hcrmiticity, and being positive semi-definite. 

Two photon Overlap — In the analysis of the previ- 
ous section it was assumed that the photons arc indis- 
tinguishable in all degrees of freedom (except for polar- 
ization). The actual experiment, however, involves the 
interference of a CS and SPDC which in general might 
not have the same spatial and spectral modes. In the 
experiment wc use a single mode fiber to ensure spa- 
tial overlap and a band-pass filter for enhanced spectral 
overlap. To verify the validity of the indistinguishability 
assumption we measure the two-photon overlap between 
the two sources (defined below), which is the lowest or- 
der of interference possible since the lowest non- vacuum 
contribution of the SPDC has two photons. Recently, 
heralded single photons from SPDC were shown [22[ to 
overlap a single photon from CS with high visibility in a 
Hong-Ou-Mandel interference [23|. An overlap of unity 
implies perfect indistinguishability between the photons 
while an overlap of implies complete distinguishability. 
We note that in the case of highly distinguishable pho- 
tons the recently introduced method of hidden differences 



[24| may be used. 

Considering the overall quantum state of a system fed 
by two separate beams, SPDC and a CS, 

I*) = \^spdc) <8> \*cs) (7) 

where the quantum states of SPDC and CS can be writ- 
ten in terms of their two-photon spectral wave-functions 
$(wi,a;2) and A(u>i)A(u>2) respectively as 

\^SPDc) = — != J J (^l<fa2§{ui,U2)\U\)\u)2) (8) 

\^cs) = -^yj JdLo 1 dLo 2 A(oj- L )A(oj 2 )\oJi)\uJ2) (9) 

where U and V are proportional to the 2-fold coincidence 
of the SPDC and CS respectively and serve as normaliza- 
tion factors. We superpose the two states on two input 
ports of a beam splitter. The 2-fold coincidence in one of 

(2) 

the output ports of the beam splitter, denoted by P bi , 
will be proportional to 



djidw 2 |$(wi, w 2 )| 2 + J J dujiduj 2 \A(uji)A(uj2)\ 2 + 

^idw 2 dwidw2(e' 2¥,cs ^(wi)A(a;2)$*K,^) + e-^ cs A*{uj 1 )A*{u) 2 )^' x ,u)' 2 )) 



= U + V + 2VUV cos(2fcs)\(^SPDc\^cs) 
I 



(10) 



where tpes is the CS phase (relative to the SPDC). 
From Eq. (fTU|) wc can see that the last term is an 
oscillating cosine function of ifes, while the other two 
terms contribute a constant background. Therefore, by 



(2) 

measuring P£ , U, and V, and scanning the phase (pes, 

(2) 

we can fit Pi to a cosine function with a constant 
1 

background and extract the overlap of the normalized 
SPDC and CS states, \{^spdc\^Cs)\, which can result 
any value between (orthogonal states) and 1 (identical 
states). 



III. EXPERIMENTAL RESULTS 

The experimental setup consists three consecutive 
parts: state preparation, state transformation, and de- 
tection. NOON states are prepared for different values 
of N by setting pair amplitude ratio between the CS and 
SPDC inputs denoted by 7, and thephase of CS relative 
to the SPDC represented by ip C s at the ±45° X, Y 
polarizations basis respectively (see FigJT]). A transfor- 
mation on the prepared states is performed using a phase 



shifting LC and a rotated HWP to (iV+1) 2 different pairs 
of angles (<fi,0), which are given in Tabic Q] The result- 
ing states are measured using five single photon counting 
modules (SPDC) in polarization mode b\. 

We have also performed a theoretical simulation of 
the experiment using the definitions of a CS and SPDC 
[l3| . The simulation accounts for all sources of loss in 
the experiment in the abstract loss model [25j using a 
single parameter 77 representing the overall transmission 
of the setup. The simulation also assumes that all the 
photons arrive in a single identical spatial and temporal 
mode. Good agreement of the simulation and the exper- 
iment implies that the assumption of indistinguishability 
is valid. 

We start by measuring the two photon overlap (TPO) 
function (defined in the previous section) which describes 
the two-photon interference of the CS with the SPDC on 
a beam splitter. To perform this, in the actual experi- 
ment, we use our setup (see Fig. [I}, while 6 = 0, <fi = tt/2 
serves as a beam splitter mixing the SPDC and CS, and 
measure the 2-fold coincidence in polarization mode 61, 
one of its output ports. The LC in the CS path con- 
trols the phase, tpes, which enables the TPO fringes. 
Fig. [5] shows the measurement of the two photon over- 
lap, for SPDC and CS 2-fold coincidence rate of about 
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FIG. 2: (Color online) Two photon overlap measurement of 
CS and SPDC. 2-fold coincidence measurement of an out- 
put port of a beam splitter, Pbm, fed by CS in one of its 
input port and SPDC in the other. Error bars indicate 
±<t statistical uncertainty. Solid line is obtained by fitting 
the data according to Eq. (|10[) resulting with an overlap of 
|(*spdc|*cs)| = 0.970±0.020 (with 95% confidence bound, 
see text). 

50Hz, which is low enough to significantly neglect the 
next contributing order, which arises from a three pho- 
ton event with one photon loss. As follows from Eq. (fTU)) 
by fitting the measurement to a cosine function of ipcs 
with a constant background (2-fold coincidence rate of 
SPDC and CS) the TPO between the CS and the SPDC 
is extracted, \(^spdc\^cs)\ = 0.970 ± 0.020 (with 95% 
confidence bound). 

The experimental reconstructed density matrices for 
the prepared states for N = 2, 3, 4 are shown in Fig. [3^. 
The fidelity parameter T of an estimated state p es t to 
overlap a theoretical state pth is defined as 




It gives a quantitative measure of the correspondence 
between theoretical and experimental density matrices of 
the states. We compare the experimental results both to 
the theoretical NOON states, and to the theoretical simu- 
lation of our setup. The experimental error in the fidelity 
is calculated from a MC simulation of the statistical error 
of the TV-fold coincidence rate and applying the density 
matrix reconstruction procedure. The fidelities of the 
density matrices of NOON state for N = 2, 3, 4 (Fig. [3ji) 
to the theoretical ideal NOON states arc 0.931 ± 0.004, 
0.895 ± 0.021 and 0.721 ± 0.063 respectively, while their 
fidelities to the states resulting from the simulation are 
0.984±0.003, 0.986±0.014 and 0.930±0.089 respectively. 
The high fidelity of the experiment to the simulation in- 
dicates that the main source of deviation from an ideal 
NOON state indeed stems from the low efficiency of our 
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FIG. 3: (Color online) Experimental results: real and imag- 
inary parts of density matrices of (a) NOON states for N = 
2, 3, 4 (b) correlated photon holes of 3,0 and 0,3 (in colored 
bars), reconstructed in the photon- number basis of two modes 
states labeled \m,n), where m and n are the number of pho- 
tons in polarization modes fii and a,2 respectively (see Fig. [1}. 
Solid frames represent density matrices of the ideal states. 



setup and not from photon distinguishability of the two 
photon state. 

Fig. [4] shows the dependence of the fidelity (of NOON 
states to ideal NOON states) on r/ (defined as the coin- 
cidence to singles ratio of the SPDC, being the overall 
setup transmission,), for N = 2, 3, 4. A simulation of 
our system (as described above) for different values of r\ 
with our procedure of density matrix reconstruction was 
employed. The low fidelity for lower transmission shows 
the susceptibility of the high-NOON states to loss. The 
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FIG. 4: (Color online) NOON states fidelity dependence on 
setup transmission r\. Fidelity of the states generated by our 
setup to ideal NOON states are shown as a function of the 
setup transmission r\ for iV = 2, 3, 4 in blue, red and green 
respectively. Simulations (solid lines) are shown together with 
the measurements (squares) in the current setup transmission, 
77 = 0.11. Error bars indicate ±u statistical uncertainty. 

fidelity in our setup is shown here to degrade mostly by 
the setup transmission. Our current setup transmission, 
r) = 0.11, is mainly limited by three major technicalities: 
detectors efficiency, coupling of light into the single-mode 
fiber, and the band-pass filter transmission. 

Additionally, we show in Fig. density matrix recon- 
struction for yet another class of quantum states of light, 
the CPH states [l2| • These states exhibit absence of mul- 
tiphoton coincidence events between two modes amid a 
constant background. We show here as an example the 
state reconstruction of a CPH with holes in the states 
|3; 0)^3 and |0, 3)^^, where A and B correspond in our 



setup (Fig. [IJ to polarization modes a\ and 0,2 respec- 
tively. This is achieved by setting 7 = 3 and ipcs — n - 
The fidelity of the measured state to the ideal CPH state 
is 0.888 ± 0.037. 



IV. CONCLUSION 

We have performed a complete analysis of High-NOON 
states generated by mixing CS with SPDC. We recon- 
structed density matrices for the measured states show- 
ing high fidelities with respect to other methods [H, [H| 
in a relatively simple setup with no background substrac- 
tion. The quantum state tomography was performed un- 
der the simplifying assumption that the photons are in- 
distinguishable in all degrees of freedom but polarization. 
We checked this assumption by measuring the two pho- 
ton overlap. The results were analyzed in terms of fidelity 
with ideal NOON states, and a simulation, and are lim- 
ited in our configuration by the loss which is dominated 
by three elements: single mode fiber, band pass filter, 
and single photon detectors, which roughly contribute 
an amount of 50% loss each. We thus stress the need in 
biphoton sources which can be spectrally and spatially 
matched to CS and single photon detectors with high 
quantum efficiencies. 
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